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Abstract. Coupling between the stratosphere and the tropo- 
sphere allows changes in stratospheric ozone abundances to 
affect tropospheric chemistry. Large-scale effects from such 
changes on chemically produced tropospheric aerosols have 
not been systematically examined in past studies. We use 
a composition-climate model to investigate potential past 
and future impacts of changes in stratospheric ozone de- 
pleting substances (ODS) on tropospheric oxidants and sul- 
fate aerosols. In most experiments, we find significant re- 
sponses in tropospheric photolysis and oxidants, with small 
but significant effects on methane radiative forcing. The re- 
sponse of sulfate aerosols is sizeable when examining the 
effect of increasing future nitrous oxide (N 2 O) emissions. 
We also find that without the regulation of chlorofluorocar- 
bons (CFCs) through the Montreal Protocol, sulfate aerosols 
could have increased by 2050 by a comparable amount to the 
decreases predicted due to relatively stringent sulfur emis- 
sions controls. The individual historical radiative forcings 
of CFCs and N 2 O through their indirect effects on methane 
(—22.6 mW m -2 for CFCs and —6.7 mW m -2 for N 2 O) and 
sulfate aerosols (—3.0 mW m -2 for CFCs and +6.5 mW m -2 
for N 2 O when considering the direct aerosol effect) dis- 
cussed here are non-negligible when compared to known 
historical ODS forcing. Our results stress the importance 
of accounting for stratosphere -troposphere, gas-aerosol and 
composition-climate interactions when investigating the ef- 
fects of changing emissions on atmospheric composition and 
climate. 


1 Introduction 

Stratospheric ozone (O 3 ) changes have been large over the 
past few decades mainly due to anthropogenic emissions of 
ozone depleting substances (ODS), and are expected to be 
large in the future due to the phase-out of ODS emissions 
and to climate change (e.g. Austin et al., 2010; Eyring et al., 
2010). Stratospheric ozone exerts direct radiative forcing on 
the Earth system (Forster et al., 2007), while it also indirectly 
affects other climate aspects such as tropospheric circulation 
and precipitation (Son et al., 2008; Kang et al., 2011). Fur- 
thermore, its changes can impact other radiatively important 
species, especially by affecting the tropospheric oxidizing 
capacity. The latter is largely determined by the abundances 
of hydroxyl radicals (OH), which depend on tropospheric 
ozone and nitrogen oxides (NO x ), humidity, abundances of 
species that consume OH (carbon monoxide (CO), methane 
and non-methane volatile organic compounds (NMVOCs)), 
and ultra-violet (UV) radiation (e.g. Thompson, 1992; Spi- 
vakovsky et al., 2000; Lelieveld et al., 2002; Voulgarakis et 
al., 2013). 

Changes in stratospheric ozone can affect tropospheric ox- 
idants via stratosphere-troposphere exchange (STE) (Collins 
et al., 2003; Hegglin and Shepherd, 2009; Prather and Hsu, 
2010; Zeng et al., 2010; Voulgarakis et al., 2011), which en- 
riches the troposphere with ozone, and by driving changes 
in the photolysis of tropospheric species (Madronich and 
Granier, 1992; Bekki et al., 1994; Fuglestved et al., 1994; 
Krol and van Weele, 1997; Isaksen et al., 2005; Lamarque et 
al., 2010a; Tang et al., 2011), which depends on UV radia- 
tion and, thus, the amount of overhead ozone (e.g. increased 
stratospheric ozone columns can lead to less tropospheric 
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OH). Photolysis has a central role in driving the formation of 
OH (Madronich and Granier, 1992; Fuglestved et ah, 1994; 
Thompson, 1992; Rohrer and Berresheim, 2006; Voulgarakis 
et ah, 2009), since it leads to the breakdown of ozone, 

0 3 + hv->- O 2 +O( 1 D)(k<320nm) (Rl) 

and the excited oxygen atoms produced react with water va- 
por to produce OH: 

0(*D) + H 2 0 — »■ OH + OH. (R2) 

Apart from cleansing the atmosphere from pollutants such as 
CO and NMVOCs, oxidation by OH also affects methane 
and sulfate aerosols, being the major sink of the former 
(Madronich and Granier, 1992; Voulgarakis et ah, 2013) and 
a major chemical source of the latter via gaseous oxida- 
tion of sulfur dioxide (SOi) (Unger et ah, 2006; Shindell et 
ah, 2009). Thus, OH plays a significant indirect role in de- 
termining the radiative balance of the troposphere. Sulfate 
aerosols are also produced through SO 2 aqueous oxidation 
by hydrogen peroxide (H 2 O 2 ) (Koch et ah, 2006; Unger et 
ah, 2006), and the production/loss of the latter also heavily 
depends on photolysis (e.g. Lelieveld et ah, 2002). Sulfate 
aqueous-phase production is globally more efficient than the 
gas-phase process (via OH), but the former produces aerosols 
with shorter lifetime, due to the usual simultaneous presence 
of clouds. 

Although effects of stratospheric ozone changes on tropo- 
spheric oxidants have been studied in the past (see above), 
their possible indirect link with tropospheric aerosols and 
their radiative forcing has not been examined. Here, we sys- 
tematically explore these links in a number of model exper- 
iments, considering both past and future ODS changes. We 
further investigate such stratosphere-troposphere interactions 
in a hypothetical future scenario where ODS have not been 
regulated by the Montreal Protocol. Past studies examining 
this topic did not investigate tropospheric effects (Morgen- 
stern et al., 2008; Newman et al., 2009). In Sect. 2, we outline 
the main characteristics of the model used, and describe the 
model experiments performed. Section 3 presents and dis- 
cusses the results of the experiments, while Sect. 4 outlines 
the conclusions of this study. 

2 Model description and experimental set-up 

2.1 Model description 

We used a version of the GISS composition-climate model 
(Shindell et al., 2006, 2013a), which is almost identical to 

the ModelE2 version developed for use in simulations in 
support of the Intergovernmental Panel on Climate Change 
(IPCC) Fifth Assessment Report (AR5). The benefit of us- 
ing a composition climate-model as opposed to a chemistry- 
transport model (CTM) is that the radiative forcing of the 
perturbations applied can be diagnosed immediately, with- 
out the need of extra radiative transfer calculations. The 


model has a 2° latitude by 2.5° longitude horizontal reso- 
lution and 40 vertical layers from the surface to 0. 1 hPa. Tro- 
pospheric/stratospheric chemistry includes 156 chemical re- 
actions among 51 gas species. Photolysis rates are simulated 
using the Fast-J2 scheme (Wild et al., 2000; Bian and Prather, 
2002), which accounts for the effects of modelled overhead 
ozone, clouds, aerosols and surface reflections. Aerosol trac- 
ers include sulfate (Koch et al., 2006, 2007), nitrate (Bauer 
et al., 2007), elemental and organic carbon, sea salt (Koch et 
al., 2006) and dust (Miller et al., 2006), while we also include 
aerosol indirect effects based on Menon et al. (2008, 2010). 
The model’s sulfate aerosol scheme includes prognostic sim- 
ulations of the mass distributions of dimethyl sulfide (DMS), 
methanesulfonic acid (MSA), SO 2 , and sulfate (Koch et al., 
2006). The sulfate lifetime in the model is equal to 3 days. 

The model’s skill in capturing tropospheric oxidants and 
sulfate has been evaluated in past studies and shown to be 
realistic (Koch et al., 2006; Shindell et al., 2006). A recent 
model evaluation, in which a very similar version of the 
model was used (Shindell et al., 2013a), showed that the 
model has improved significantly compared to its previous 
version (Shindell et al., 2006) when it comes to ozone sea- 
sonality in the stratosphere. However, the Antarctic ozone 
hole extends more equatorward than seen in the observations, 
and it lasts somewhat too long. Still, the total ozone biases 
in middle and low latitudes, which are of particular interest 
in this study (due to the abundance of oxidants and aerosols 
in these regions), are small, and typically less than 5 %. In 
addition, the model’s aerosol performance has thoroughly 
been evaluated against surface and satellite observations as 
part of the Atmospheric Chemistry and Climate Modeling 
Intercomparison Project (ACCMIP; Shindell et al., 2013b). 
The model’s performance was fairly satisfactory, with bi- 
ases against AeroNet measurements being among the lowest 
in the multi-model set, particularly over Europe and North 
America. This is important, as these regions are key for the 
responses we are examining. Performance was less good e.g. 
over North Africa (overestimates) and over southern midlati- 
tude regions (underestimates), however such regions are less 
crucial for the current analysis, as sulfate aerosols are not 
abundant in them. 

The global tropospheric mean mass-weighted OH concen- 
tration is 9.9 x 10 2 * * 5 molecules cm -3 , which is in good agree- 
ment with past estimates (Lawrence et al. (2001) found a 
range of 8.8-11.4 x 10 s molecules cm -3 ). Furthermore, it 
has been shown (Voulgarakis et al., 2013) that the model’s 
tropospheric chemical methane lifetime (10.6 yr) agrees 
very well with a recent observationally constrained estimate 
(11.2 ± 1.3 yr; Prather et al., 2012). The above indicate that 
tropospheric oxidising capacity is simulated reasonably in 
the model. The model’s sulfate has been found to be more 
sensitive to OH changes than several other models (Fry et al., 
2012), though these other models did not include liquid tracer 
budgets (i.e. gas phase and dissolved tracers are not tracked 
separately), which may mean that they are missing some of 
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the physics that are included in the GISS simulations, though 
differences in sulfate chemistry could be important as well. 

2.2 Experimental set-up 

In the experiments presented here, present-day anthro- 
pogenic and biomass burning emissions of NO x , CO, 
NMVOCs and aerosols follow Lamarque et al. (2010b), 
which were created in support of the Coupled Model In- 
tercomparison Project Phase 5 (CMIP5). Lightning NO x 
emissions depend on the model’s convection (Price et al., 
1997; annual present-day mean: 6.9 TgNyr -1 ), while iso- 
prene emissions respond to temperature and solar radi- 
ation based on the algorithm of Guenther et al. (2006) 
(annual present-day mean: 467Tgyr -1 ). Natural (non- 
agriculture) emissions of NO x from soils are prescribed at 
fixed values (2.7TgNyr -1 , as in Shindell et al., 2013a; 
total soil emissions including the agricultural source are 
equal to 3.4 TgNyr -1 ), as are biogenic emissions of 
alkenes (16TgCyr -1 ), paraffins (14TgCyr -1 ), and ter- 
penes (192 Tgyr -1 ). The concentrations of long-lived gases, 
such as CO 2 , methane, N 2 O and CFCs, are prescribed at the 
surface to observed present-day values. When perturbed in 
our experiments (CFCs, N 2 O; see next paragraph), they ei- 
ther follow past observations, or future projections according 
to the Representative Concentration Pathway (RCP; van Vu- 
uren et al., 2011). The RCP projections for N 2 O and CFC 
concentrations are based on calculations with the reduced- 
complexity coupled carbon cycle climate model MAGICC 

6.3 to estimate mixing ratios that would result from the RCP 
emissions (Meinshausen et al., 2011). 

We performed six main simulations (see Table 1), in all 
of which climatic conditions were kept to present-day lev- 
els: (a) BASE, which represents non-perturbed present-day 
(2000) conditions; (b) 1960CFCs, in which CFCs are set 
to 1960 levels; (c) 2100CFCs, in which CFCs are set to 
2100 levels following the RCP8.5 scenario; (d) WACFCs, 
in which CFCs are set to 2050 levels following a “world 
avoided” scenario, where CFCs are not regulated by the 
Montreal Protocol (Newman et al., 2009); (e) 1850N2O, 
in which N 2 O is set to 1850 levels; and (f) 2100N20, in 
which N 2 O is set to 2100 levels (following RCP8.5). Non- 
ODS emissions/concentrations and sea-surface temperatures 
(SSTs) are representative of present-day in all experiments 
in order to isolate the ODS effects alone. Note that allow- 
ing ODSs to affect SSTs could have a significant impact 
on the results, since both CFCs and N 2 O are important cli- 
mate forcers. In addition to the above six experiments, we 
also performed a present-day simulation with 2050 RCP4.5 
aerosol emissions, for comparison. In each experiment, we 
ran the model using the same conditions for 7 yr and then 
averaged in order to increase the signal-to-noise ratio. Two 
years of spin-up were performed in advance, which we con- 
firmed were enough for the stratospheric concentrations to 
fully adjust to the perturbations applied at the surface. 


Table 1 . Percentage (%) differences in surface CFC and N 2 O con- 
centrations between the various sensitivity experiments and the 
BASE simulation. 



CFCs 

n 2 o 

BASE 

0 

0 

1960CFCs 

-94 

0 

2100CFCS 

-76 

0 

WACFCs 

+ 1026 

0 

1850N2O 

0 

-13 

2100N20 

0 

+38 


3 Results and discussion 

Figure 1 shows the annual zonal mean photolysis rate of 
ozone to yield a singlet oxygen atom (.1(0 1 D); see Reaction 
Rl), as well as its changes in the various sensitivity experi- 
ments. .1(0 1 D) is largest in the tropics, where solar zenith an- 
gles and overhead ozone columns are smallest. We find that 
all the perturbations applied to ODS have sizeable effects on 
tropospheric photolysis (Fig. 1). The changes in .1(0 1 D) are 
largest over polar regions, where stratospheric processes de- 
pleting ozone are more effective, and thus UV changes af- 
fecting tropospheric photolysis are also large. However, even 
in non-polar regions, such as northern midlatitudes, there 
are sizeable changes (mostly larger than 5 %), which can 
have important impacts on tropospheric chemistry. The ef- 
fects are dramatic in the WACFCs experiment, where mid- 
latitude .1(0 1 D) increases by more than 120 %. 

Figure 2 shows differences in total ozone column, surface 
J(0 1 D), surface OH and surface sulfate concentrations, be- 
tween the WACFCs and the BASE simulation, for the north- 
ern extratropics, the region where sulfate is most important 
for climate and air quality. Total ozone for 2050 is lower by 
50-60 % due to the lack of regulation in CFCs in this sce- 
nario, in fairly good agreement with previous work (New- 
man et al., 2009). This leads to ubiquitous surface J(O'D) 
increases of more than 150% in the northern extratropics, 
which, in turn leads to OH increases of 40-80 %. The high- 
est OH increases are located in the regions with the largest 
total ozone and 1(0*0) changes combined with high OH re- 
cycling efficiencies due to high amounts of NO x (Eastern US, 
Western Europe, East Asia) (Lelieveld et al., 2002). 

We find subsequent large increases of sulfate aerosol con- 
centrations throughout the northern extratropics (Fig. 2d). 
Differences range from 20 to 50 % in highly populated mid- 
latitude regions where aerosols cause serious human health 
damage. The increases are around 30-35% over the north- 
ernmost states of the US, over 40 % over the British Isles, 
and 15-30 % over East Asia. Increases in sulfate aerosols are 
correlated in several cases to the enhancement of the tropo- 
spheric oxidizing capacity, as reflected by OH (e.g. over Pa- 
cific and Atlantic pollution outflow regions), which leads to 
more SO 2 oxidised to sulfate. Figure 3 shows the zonal mean 
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latitude latitude 



latitude latitude 

Fig. 1. Photolysis rate changes in the sensitivity experiments. Annual zonal mean J(0*D) in the BASE simulation (a), and difference of 
annual zonal mean J(0*D) in the 1960CFCs (b), the 2100CFCs (c), the WACFCs (d), the 1850N2O (e) and the 2100N20 (f) simulations, 
from the BASE run. Note the difference in the colour scale in (d). 


differences in concentration and chemical gas-phase produc- 
tion rate of sulfate, between WACFCs and BASE. It is clear 
that areas with large chemical production changes are also 
areas with high concentration changes, the “hot-spot” being 
the northern midlatitudes at low altitudes. This finding is in 
qualitative agreement with an earlier study on a related topic 
(Toumi et al., 1994). 

Note that changes in H 2 O 2 and aqueous-phase production 
of sulfate are also detectable in our perturbation experiments, 
especially in WACFCs, where global mean H 2 O 2 changes 
by around 38%, while global mean sulfate aerosol produc- 
tion via H 2 O 2 changes by 1 1 % (compared to 28 % for gas- 
phase), relative to BASE. However, as mentioned earlier, the 
gas-phase production changes (via OH) are expected to be 
more important due to the longer lifetime of the sulfate pro- 
duced via gaseous processes. 

Table 2 summarizes the effects of the perturbations ap- 
plied in the different experiments on tropospheric OH and 


surface sulfate concentration, as well as on radiative forc- 
ing by methane and sulfate. Note that sulfate radiative forc- 
ing calculations only take the direct aerosol effect into ac- 
count. We focus on global effects, and on effects in north- 
ern midlatitudes, where sulfate is important, both for air 
quality and climate. The OH differences are significant in 
all the experiments, with the WACFCs scenario showing 
clearly the largest effects (more than an order of magnitude 
larger than in other simulations), consistent with photolysis 
changes (Fig. 1). Global methane radiative forcing strongly 
decreases in WACFCs (by ~ 0.5 W m -2 ) due to tropospheric 
OH increases, while there are also significant changes in the 
other experiments (6-23 mW m -2 ). 

For sulfate, in terms of both concentration and forcing, 
the two cases for which we find the most sizeable relative 
changes are WACFCs and 2100N20. In the former, the dif- 
ferences range from around 12 to 20%, while in the lat- 
ter, they are around 1.7-2. 5%. Regionally, concentration 
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Fig. 2. Differences in annual mean total ozone column (a), J(0*D) 
(b), OH (c) and sulfate concentration (d), between WACFCs and 
BASE. All differences are in percent (%), and for J(0*D), OH and 
sulfate they refer to the surface. Solid lines in (a) denote the midlat- 
itude region (20-60° N) analysed in Table 2. 


differences can be as large as 50% in WACFCs and 15 % in 
2100N20, and forcing differences reach 35% in WACFCs 
and 13 % in 2100N20 (see Fig. 2d and Fig. 4), with a fairly 
good correspondence of high response regions in the two sce- 
narios (e.g. north Atlantic, East Asia, northwestern Pacific). 
Areas affected heavily include East Asia, parts of Europe and 
parts of Northeastern America and the North Atlantic. Note 
though that for 2100N20, the radiative forcing differences 
(Fig. 4a) are only significant where the signal is relatively 
strong (e.g. East Asia). An unexpected feature is that for the 
1850N2O simulation, sulfate is higher than in BASE (see Ta- 
ble 2), despite the fact that there is less OH in the atmosphere. 
However, this effect is rather small and is likely related to 
meteorological feedbacks resulting from atmospheric forc- 
ing, or to unforced variability. 

In an additional experiment where sulfur emissions were 
set to 2050 values following the R.CP4.5 scenario (in which 
SO 2 emissions decrease by more than 50 %) but ODSs were 
kept to present-day levels, we found global surface sulfate 
aerosol concentration reductions of 24 %. This indicates that 
the indirect effect of ODS on aerosols discussed here (17% 
globally) could have cancelled out part of the benefit gained 




-60 -40 -20 0 20 40 60 

latitude 


Fig. 3. Annual zonal mean differences in sulfate concentration (a) 
and gas-phase production (via OH) (b), between WACFCs and 
BASE. 


from aerosol emissions reductions, had the Montreal Pro- 
tocol not been implemented. Another interesting feature is 
that several of the changes that we document here would 
likely cancel out, at least partially, in a future scenario in- 
cluding both CFC and NoO changes (compare 2100CFCs 
and 2100N20 columns in Table 2). However, in such a case 
the dominant driver of future changes would be expected to 
be the large sulfur emission decreases in scenarios such as 
RCP8.5. Note that total sulfate burden changes depend on 
the current relative amount of SO 2 removed by deposition 
versus conversion to sulfate, thus it would be useful to test 
the variation of this balance across different models in fu- 
ture work (in our model, the SO 2 deposited globally is about 
double the amount that is chemically converted). 

The indirect radiative forcing of ODSs through the pro- 
cesses presented here has not been discussed in previous 
studies (Forster et al., 2007). We compared these effects to 
the magnitude of previously estimated direct historical ODS 
forcing, as well as to indirect forcing through other pro- 
cesses. Our calculated global methane radiative forcing due 
to past CFC changes is equal to about 7% of the direct 
preindustrial (PI) to present-day (PD) forcing from CFCs 
(0.32 W m -2 ), and to 65 % of their indirect effect through 
stratospheric ozone radiative forcing (— 0.04 Wm -2 ), based 
on Forster et al. (2007). For sulfate radiative forcing, the 
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Table 2. Differences in annual mean tropospheric OH concentration, surface sulfate aerosol concentration, methane and sulfate aerosol 
instantaneous top-of-the-atmosphere (TOA) radiative forcing (RF), between the sensitivity experiments and the BASE run 3 . 



1960CFCs 

2100CFCs 

WACFCs 

1850N2O 

2100N20 

GL b trop. 
OH conc. c 
(%) 

-3.0 ±0.2 

-2.7 ±0.2 

±130.9±0.8 

-0.9 ±0.2 

±2.2 ±0.2 

OL 

CH 4 RF d 
(mW m -2 ) 

±22.6 ± 1.0 

±20.4 ± 1.5 

-515.1 ± 1.4 

±6.7 ± 1.1 

—16.2 ±1.0 

NM fo trop. 
OH cone. (%) 

-2.1 ±0.1 

— 1.8 ±0.2 

±80.5 ± 1.0 

— 1.0 ±0.2 

±2.7 ±0.2 

GL surf. 

Sulfate conc. c (%) 

-0.3 ±0.3 

-0.8 ±0.3 

±17.5 ±0.3 

±0.5 ±0.4 

±1.7 ±0.4 

NM surf. 

Sulfate cone. (%) 

±0.4 ±0.4 

-0.5 ±0.3 

±19.9±0.4 

±1.1 ±0.5 

±2.4 ±0.6 

GL 

Sulfate RF c e 
(mWm -2 ) 

±3.0± 1.6 

±6.0±2.2 

—83.9 ±1.7 

-6.5 ±1.6 

-12.1 ±2.6 

NM 

Sulfate RF 
(mW m -2 ) 

-3.7 ±2.8 

±2.6±6.5 

—157.2 ±4.1 

— 17.4±6.0 

— 18.3 ± 7.3 


a Seven years from each simulation were averaged in order to increase the signal-to-noise ratio. The uncertainties are ±l<r standard 
deviations of the mean, based on the fluctuations among the different years in each simulation. CFCs and N 2 O concentrations for 
2100 follow the RCP8.5 scenario. 
b GL = global; NM = northern midlatitudes (20-60° N). 

c Mean tropospheric OH concentrations used for calculating the differences are mass-weighted, while mean surface sulfate 
concentrations and mean sulfate radiative forcings are area-weighted. 

d We infer approximate global methane changes offline to calculate methane radiative forcing changes presented in Table 2. Global 
methane is calculated from global OH changes assuming a liner relationship between OH and methane chemical loss, and an 
additional chemical feedback (<$ln(OH)/<$ln(CH4) = —0.32; Prather et al., 2001). Methane radiative forcing was calculated 
following the standard formulations used in the IPCC Third Assessment Report (TAR; Prather et al., 2001). 
e Sulfate direct instantaneous radiative forcing is calculated online as average over time during the simulations, with radiative 
transfer calculations performed twice at each point in time with the only difference in the two calculations being the sulfate field. 


CFC effects here are smaller and less significant (1 % of 
the direct CFC forcing and 8% of the their indirect forc- 
ing through stratospheric ozone changes). The historical in- 
direct effects of N 2 O on radiative forcing through methane 
and sulfate (roughly equal) are equal to around 5 % of its di- 
rect PI-PD forcing (0.15 W m" 2 ) and to 70 % of its indirect 
effects through stratospheric ozone changes (—0.01 W m~ 2 ). 
The effects of future N 2 O changes on radiative forcing in 
our experiments are 2-3 times larger than the historical ef- 
fects (due to the larger increases in 2000-2100 compared to 
1850-2000), while the future effects of CFCs are also equal 
or larger to the historical. 

Changes in stratospheric ozone can impact tropospheric 
oxidation via STE, and this effect is expected to be oppo- 
site of the photolysis effect. Less ozone in the stratosphere 
can lead to less ozone transported to the troposphere and, 
thus, less OFI produced by ozone photolysis. Flowever, the 
effect of UV changes on photolysis appears to be dominat- 
ing the OF! responses that we find in the troposphere. There 
is even one case (WACFCs, Southern Flemisphere) in which 
ozone depletion in the stratosphere is followed by ozone in- 
creases in the troposphere. Flowever, this is driven by large 
increases of tropical/ southern sub-tropical stratospheric NO x 


a 


Q 

=! 

05 



Change in sulfate aerosol shortwave radiative forcing (2100N20 - BASE) 



Fig. 4. Changes in sulfate radiative forcing. Annual mean difference 
of midlatitude sulfate aerosol radiative forcing (shortwave, top-of- 
the-atmosphere) in the 2100N20 (a) and the WACFCs (b) simula- 
tions, from the BASE run. Note that changes are statistically signif- 
icant only for isolated regions in the mid-north Atlantic, Mediter- 
ranean and northern East Asia in (a), but for almost all mid-latitude 
regions in (b). 
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subsequently transported into the Southern Hemisphere tro- 
posphere (also see Fig. Id for effect of this feature on photol- 
ysis). It remains unclear whether this feature is realistic, as 
stratosphere-troposphere chemical interactions in a “world- 
avoided” scenario have not been examined in other stud- 
ies. Note that we find that tropospheric temperature and hu- 
midity changes are very minimal (< 1 %) in all the experi- 
ments (which is expected since we have fixed SSTs/sea ice 
to present-day values), so we rule out the possibility that 
climate changes induced by differences in CFCs and strato- 
spheric ozone drive any of the changes simulated. 

Our experiments presented here are idealised, since a sce- 
nario in which only ODS will change, while other aspects 
of the climate system remain constant is not realistic. Nev- 
ertheless, by isolating the ODS effect, we provide insight 
into particular processes which can be important. For the 
two scenarios in which we found the most sizeable differ- 
ences (WACFCs and 2100N20), changing future emissions 
of gases that affect tropospheric oxidation chemistry and cli- 
mate change would be expected to add to the effects of ODS 
on aerosols. In particular, studies have found that continen- 
tal regions, where sulfate aerosols are of high importance, 
will likely experience increases in OH abundances due to in- 
creases in NO x emissions (e.g. SRES A1B and A2 scenar- 
ios; see Nakicenovic et ah, 2000), while the global temper- 
ature and humidity increases expected due to global warm- 
ing may lead to generally more OH production in the tro- 
posphere (Wild and Palmer, 2008). This, in addition to the 
oxidant increases due to ODS effects found here, could lead 
to even larger sulfate responses in WACFCs and 2100N20. 
Note though, that when more optimistic scenarios are used 
for future NO x emissions (such as in the RCP scenarios; van 
Vuuren et ah, 2011), there is a mixture of positive and neg- 
ative OH responses over industrialised regions, with large 
inter-model variability (Voulgarakis et ah, 2013). Further- 
more, both sulfur and NO x emissions may shift closer to the 
tropics in the future, which implies that some of the effects 
at lower latitudes may be more important than we find here. 

Projected increases in tropical upwelling could lead to 
further lower stratospheric ozone decreases at low latitudes 
(Eyring et ah, 2010), which would, again, add to the UV, ox- 
idant and sulfate aerosol responses. Also, projected changes 
in atmospheric circulation may drive increased transport of 
ozone from the stratosphere to the troposphere (Collins et 
ah, 2003; Shindell et ah, 2006), increasing tropospheric OH 
production. However, in a scenario like 2100N20, if we ac- 
count for the simultaneous decline in CFCs, part of the tro- 
pospheric effects will be cancelled out. Finally, note that in 
the calculations of OH and sulfate responses, feedbacks on 
OH through changes in methane have not been taken into 
account, since methane concentrations were kept constant. 
Had these feedbacks been included, all responses presented 
here would have been stronger. For example, only due to the 
induced methane decreases in WACFCs, global OH would 
have dropped by 36 %. 


4 Conclusions 

We performed model experiments to investigate the role 
of ozone depleting substances in driving changes in tropo- 
spheric oxidants, aerosols and radiative forcing. These links 
have not been examined in previous studies investigating the 
effects of ODS or the evolution of tropospheric aerosols. Our 
results indicate that changes in ODS and stratospheric ozone 
can induce large-scale changes in tropospheric composition, 
with a non-negligible effect on radiative forcing. Indirect his- 
torical forcings of CFCs and N 2 O through their effects on 
tropospheric methane and sulfate range from 8% to 70% 
(depending on the species) of the known corresponding in- 
direct forcings through stratospheric ozone depletion. Spe- 
cial attention is needed regarding NtO effects beyond 2100, 
since its emissions are expected to be the dominant driver 
of future stratospheric ozone depletion (Ravishankara et ah, 
2009). However, we also note that the future effects of N 2 O 
and CFCs are of opposite signs and may thus partly cancel 
out. 

Our analysis also suggests that CFCs emissions reduc- 
tions resulting from the Montreal Protocol may be trans- 
lating to tropospheric composition and climate changes that 
have not been discussed in assessments exploring the “side- 
effects” of the agreement (Velders et ah, 2007). We suggest 
that the Montreal Protocol may have been a factor mediat- 
ing sulfate aerosol concentrations, especially when consider- 
ing effects later in the 21st century. Finally, our results im- 
ply that accounting for stratosphere-troposphere, gas-aerosol 
and composition-climate interactions (not fully considered 
here) is important when assessing the impact of regulations 
in emissions. Earth system models, which are designed to in- 
clude all the above-mentioned interactions, are expected to 
facilitate such efforts in the future. 
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